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Abstract 

The role of double humped states in spreading of wave packets for the nonlinear Schrodinger 
Q\ equation (NLSE) with a random potential is explored and the spreading mechanism is unraveled. 
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Comparison with an NLSE with a double-well potential is made. There are two independent 

(N 

affects of the nonlinearity on the double humped states for the NLSE: coupling to other states and 
^ destruction. The interplay between these effects is discussed. 
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We consider the discrete nonlinear Schrodinger equation with a random potential in one 
dimension: 

= -4>n+l ~ ?Pn-l + tnlpn + (3 \4>n\ 2 Ipn (l) 

Where e n are random potentials chosen uniformly from the interval [—2, 2] and (3 is a pos- 
itive constant. For (3 = this is the Anderson model, where all the states are localized. 
Consequently, a wave packet that is initially localized will remain localized in the vicinity 
of its initial position. A question that is subject to extensive research is whether Anderson 
localization can survive the nonlinear term j3 \tp\ 2 ip [T]. Numerical simulations indicate that 
for Q Anderson localization is destroyed and subdiffusion takes place [21 El IH El E]- Heuris- 
tic arguments were developed in order to explain these results [H El E], but the detailed 
mechanism of possible spreading is not clear. Resonances between eigenstates of the linear 
model, namely Q with /3 = provide a reasonable mechanism for spreading and it is the 
subject of the present paper. 

Double humped states y?+,y?- are two eigenstates of the Hamiltonian ([!]) which are lo- 
calized over the same two sites that are far in real space while their energies are very close. 
An example for such states appears in Fig.l. According to Rabi's formula, in the linear case 
(f3 = 0), if one places (at time t = 0) a wave packet on the site of one hump and the states 
are exactly symmetric or antisymmetric with respect to the interchange of the humps, one 
finds the packet on the other site in time t with the probability 

P l2 (t) = sin 2 Ue~) (2) 

where AE is the difference between the energies of the two double humped states and the 
time period is T^ a u = ^j. The period is preserved also for the case when the symmetry of 
the hump interchange is broken, as in the case of the random potential. This mechanism 
of jumping between sites has proved to be the main mechanism for low frequency ac con- 
ductivity in disordered media [TJ. It is expected that its behavior may be strongly affected 
by the nonlinear term. For a double-well potential the low energy states are symmetric and 
antisymmetric double humped states, with the humps in the centers of the wells. In the 
absence of nonlinearity, ^ holds. But, for sufficiently strong nonlinearity the wavepacket 
will be confined to the initial well [HE]. In the present work we would like to explore if 
double humped states contribute to a mechanism of resonant spreading in the NLSE. 
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Figure 1: A pair of double humped states for the linear system (/? = 0). The states are marked 
with blue solid line and green dashed line. 

In our study, we would like to distinguish between resonant spreading (caused by the 
double humped states) and diffusive spreading from one state to its neighbors. For this 
purpose, we had to find realizations where the humped states are located far from each 
other (in comparison with the localization length, in our case £ ~ 6). The probability to 
find double humped states with humps located at a distance L is proportional to exp (— L/£) 
[7] and therefore realizations which couple states located far from each other are very rare (as 
we are looking for them in a finite region in real space). In order to overcome this problem and 
create a pool of realizations with double humped states in some region in real space, we have 
developed a strategy for "double humps hunting". We choose some random potential having 
localized eigenstates in the linear case (Anderson localization). We focus on two sites so that 
we will have double humped states which are localized on these two sites. These sites will 
be denoted by O and P in what follows. The Hamiltonian of this realization is diagonalized, 
which results in a diagonal matrix with eigenenergies on the diagonal. Now, we vary the site 
energy of the original model on one site (say P) of the two sites mentioned above. According 
to Feynman-Hellman theorem, when we increase monotonically the potential of a site, its 
energy is monotonically increasing and we can easily find a point where the two diagonal 
terms are approximately equal. Since we change the realization, the Hamiltonian is not 
diagonal anymore and the sites are coupled by matrix elements of the order of exp (— L/£). 
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Taking the potential realization which creates almost identical energies in the diagonal of 
the Hamiltonian (written in the initial eigenstates basis), we can usually construct double 
humped states for this realization. In this way, we found a set of realizations having double 
humped states with distance of 25 sites (about 4 localization lengths in our case) between 
the humped sites. 

After choosing appropriate realizations, we had to know which values of (3 should be 
chosen in order to see the influence of the double humps. If we choose very small values of 
(3, the system will behave similarly to the linear case and a wave packet initially localized 
on one humped site will oscillate between the humped sites for very long times. However, 
for large values of (3, the linear eigenstates become irrelevant very quickly (compared to the 
period of the oscillations) and the correlation between the double humped states is broken 
before they have a chance to affect the dynamics. Moreover, high values of (3 suppress 
the oscillations between the humped states even in the double-well case |9] where there is 
no mixing with other states. So, we have to choose the (3 values very carefully. For this 
purpose, we use a double-well model [10] to set the scale of the effect of (3. In particular 
we find numerically for each disorder realization a value of (3i for which only \ of the wave 
function oscillates between the humped sites O and P when the double humped states are 
detached in the computation from all other states. When we run the dynamics of Q for 
double humped realizations with (3 — we can see clearly the influence of the resonance 

4 

and we are still able to observe spreading for reasonable times. 

In other words, the nonlinearity has two effects: destroying the double humped states 
and populating other states of the linear model. In order to distinguish the two effects we 
compare to the double-well model with a nonlinear term, where the two lowest energy states 
can be assumed isolated from the other states. 

The difference between the double-well model and (jlj is that in the double-well model 
only two states participate in the dynamics (see appendix) and only these were taken into 
account for this model. Therefore, numerical calculations for the double-well model are 
much faster and allow us to estimate the behavior of ([I]) without performing time consuming 
(split-step) calculations. In addition, the dynamics in the double-well problem is periodic 
and gives us the time scale of the oscillations. Deviations of Q from the double-well model 
appear when additional states become involved in the dynamics. This happens, naturally, 
when we increase (3. So, first we should calculate (3 C , the largest (3 for which the double-well 
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model dynamics is still similar to Q and make sure that f3 c > f3i (otherwise, our results 
for (3i will have no clear meaning for the NLSE). We have located an initial wave packet y~o 
around one of the humps (as a superposition of the two double humped states) at site O. 
We have followed the population difference between the double humped states in the double- 
well model and in the NLSE during one time period T which is numerically calculated for 
each realization based on the nonlinear double- well model (see appendix). f3 c was defined 
as the highest (3 value for which ^ J Q T (wd ou bie-weii — wnlse) 2 dt < 0.001 where w denotes 
the population difference. [3 C is expected to be high when the overlap between the double 
humped states and other states in the system is small and we can see a correlation between 
j3 c and the parameter 

«-' = £' 



V 000 

(3) 



E — Ei 

The index i in the sum runs over all the eigenstates on the lattice except for the two double 
humped states. E is the energy of the initial wave localized at site O and Ei are the 
eigenvalues of the system. The numerator is Vf 00 = J2 n y Q n ■ v^ n where the eigenfunctions 
with center of localization at site i are denoted by vt and v ijn is the n component of vector 
while yo, n is the n component of yo. The reasoning for the importance of ^ is explained in 
[21 El E] (where R is defined in a slightly different way) and the correlation to f3 c is shown in 
Fig. 2. The correlation deteriorates when Vf 00 is replaced by other quartic combinations 
of components of y~o and T^. 

In order to see the influence of the double humped states in specific realizations, we 
should compare them to realizations where such states are broken that will be named "broken 
realizations". The "broken realizations" are the same set of realizations as the realizations 
where double humped states are found except for one fundamental difference - we have 
changed the disorder potential in one of the humped sites to be zero, namely ep = 0, and by 
this broke the coupling and destroyed the double humped states without causing qualitative 
changes to the other eigenstates of the system. For each pair of double humped and "broken" 
realizations, we have chosen an initial wavefunction y~o located around O as a superposition 
of the two double humped states and followed the evolution of the wavefunction. The 
evolution of the wavepacket in time was calculated according to Q with j3 = fti using the 
split step method [5]. A quantity which interests us when we measure the spreading of a 
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Figure 2: Correlation between R and j3 c (blue asterisks), 
wavefunction is the second moment, denned as 

m 2 = ^ (n - n) 2 \4> n \ 2 (4) 

n 

where n = |^n| 2 is the averaged location of the wavefunction. When we compare the 

growth in the second moment for double humped realizations and the broken realizations, 
we see that the second moment of the double humped realizations grows faster, when the 
realizations are selected as was outlined above and in both cases the initial wave packet 
is localized at O. Some examples are presented in Fig. 3. We examined 25 realizations 
of this form and the behavior presented in Fig. 3 is representative of all of them. This 
indicates that double humped states do substantially contribute to the spreading process of 
a wavefunction more than typical states. 

In conclusion, we see that in the presence of nonlinearity that is not too strong, the 
spreading of a wave packet prepared initially near some site O is substantially stronger if 
there is a double humped state with one of its humps near O, than if the states peaked 
near O are single humped. We found that there is a regime of values where (3 is sufficiently 
small so that the double humped structure is preserved but the packet is not only oscillating 
between the humps but also leaks to other states, leading to spreading. In order to find this 
nonlinearity regime, we have used the double-well model to isolate the two double humped 
states from the other eigenstates of (pi) with (3 — 0. We found that if (3 is small enough so 
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that the oscillations between the two states are not suppressed in the double-well model, 
then the double humped states will contribute to the spreading for the NLSE. Since double 
humped states are suppressed and do not contribute to the spreading for high nonlinearities, 
we can not conclude that they dominate the spreading of the NLSE. Exploring what is the 
dominant mechanism for this problem is left for future a research. 
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Appendix- Double-well model 

In order to predict the response of the double humped states to variations of (3, we first 
investigate a model where only two states exist, the double-well model. In this way we avoid 
the influence of the other states of the NLSE. For this model, the NLSE is 

i 9 *^ = - V 2 ^ (r, t) + [e(r)+/3|* (r, t) | 2 ] V (r, t) (5) 

where e (r) is the double-well potential. It is convenient to write the wavefunction in the 
form [TTj 

* (r, t) = V>i (t) 0i (r) + ^2 (t) 02 (r) (6) 

where 0i (r) and 02 (r) are symmetric and antisymmetric combinations of the double humped 
eigenstates (and therefore they are orthogonal) while ipi(t), -02 (£) are the amplitudes of 
0i ( r ), 02 0") at time t. Eq. ([5| takes the form: 

= - [Vi V 2 0x + ^ 2 V 2 2 ] + [e (r) + /3 |tf (r, t) | 2 ] * (r, t) (7) 
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Figure 3: The second moment as function of time for a representative double humped (solid blue) 
and broken (dashed green) realizations for wave packets started in the vicinity of O. 

After multiplying both sides by 4>i ( r ) (0i an d 02 are localized and therefore they can be 
chosen to be real) and integrating over r, Q becomes 

.dipi 



df j [^i0iV 2 0i + ^20iV 2 2 ] dr + J e [0 2 ^i + 0201^] dr 

+ (3 j dr [iPl^t + tybl>l + 2 |Vi I 2 ^2) 0?0 2 + 

+ (2 |^2| 2 ^l + V^l) 0?02 + l^2| 2 ^20^l] 



Following p2], it is convenient to write @ as 

^1 = -i (ui + |^i| 2 ) Vi -iKip 2 

-l (2A 1 ij 2 + Ay$$\ + S^'Vl + M |V>2 I 2 V>2 - 2Ai |^2 I 2 V2) 

where 



- I (|V0i| 2 + e0 2 + 2/?0 2 2 )c/r 
Sl x = -P j (4 - 2cf>l<f>l) dr 



(8) 
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K = - J (V0iV0 2 + e0!0 2 ) 
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%dr 



and we have used the relation \ipi\ 2 + \ip 2 f — 1- In a similar way, 

-% [u 2 + Q 2 l^l 2 ) ^2 - ii^V'i (11) 



(#2 



-i (2A 2 ^i + A 2 ^i + B^r 2 + A, |Vi| 2 ^i - 2A 2 |^i| 2 ^i) 



where 



^2 = - J (|V0 2 | 2 + £02 + 2/?0?02) ^ (12) 

and 

Q 2 = -pj {4 - 20 2 2 ) 

In order to establish the connection with the double humped states of Q, the coefficients 
of (10) and (12) were taken from the the Schrodinger Eq. Q. First we have expressed 



these coefficients for the linear case /3 = where only ui, u 2 and K do not vanish. For this 
purpose we find tp + and <£>_, the double humped eigenstates of ([I]) for (3 = 0. the amplitudes 
ijji (t) and i/j 2 (t) of the symmetric and antisymmetric combinations 

01 (r) = + (13) 

02 (r) = (^+ - (14) 

satisfy the Schrodinger Eqs. (|9]) and ((TT|. Therefore, When we write the Hamiltonian (jl]) in 
a basis composed from <fii and 4> 2 in addition to all the single humped eigenstates of ([!]), K 
will appear as an off diagonal term which couples 0i and (fi 2 while uj\ and cj 2 will appear as 
diagonal terms. In the nonlinear case, K stays the same while iv lj2 = u l ^ ear — 2(3 J <p\<p\dr. 
In order to find the corrections to co>i 2 for the nonlinear Hamiltonian and to calculate all the 
other coefficients ( jTo| and (12), we find numerically the vectors <pi and <p 2 with the help of 



(13) and (14) from ([I]) as explained above. 
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It is convenient to follow the dynamics described by the variables u = ipiip 2 + ^2^1, 
v = -i (-01-02 _ V^i) and w = 

IV^iI - \H ■ After some simple procedures, we obtain a 
vector equation of the motion [10], 

f^x? , 1B ) 

in which ~~p = (u,v,w) is a vector characterizing the state of the coupled system on the unit 
sphere, i.e., u 2 + v 2 + w 2 = 1, and, T = (Ti,T 2 ,T 3 ) where 

Ti = u 1 — u 2 + ^1 (1 + w) - -Q 2 (l-w) + {A x - A 2 ) u (16) 

T 2 = -Bv (17) 

T 3 = 2K + 2 (Ax + A 2 ) + fin - Ai (1 - w) - A 2 (1 + w) (18) 

It is easy to find (numerically) w (t) which gives us the time period and the amplitude of 
the double-well oscillations. For small nonlinearities, these results are good estimations of 



the NLSE behavior on a lattice. In this work, we have used (15) to find (3i values. For 

■ — > 4 

this purpose, we have chosen an initial value for w which represents a wavepacket localized 
around site O. Following the dynamics of w(t) for different values of nonlinearity /3, we have 
found the maximal (3 for which at least \ of the wavepacket is oscillating between sites O 
and P. This (3 is (3x. 
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